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ABSTRACT 



Context. The Monoceros Loop (SNR G205. 5+0.5) is a large shell-type supernova remnant located in the Rosette Complex region. It 
was suggested to be interacting with the Rosette Nebula. 

Aims. We aim to re-examine the radio spectral index and its spatial variation over the Monoceros SNR, and study its properties of 
evolution within the complex interstellar medium. 

Methods. We extracted radio continuum data for the Monoceros complex region from the Effelsberg 21 cm and 11 cm surveys and 
the Urumqi 6 cm polarization survey. We used the new Arecibo GALFA-Hl survey data with much higher resolution and sensitivity 
than that previously available to identify the HI shell related with the SNR. Multi-wavelengths data are included to investigate the 
properties of the SNR. 

Results. The spectral index a (Sy oc y") averaged over the SNR is -0.41±0.16. The TT-plots and the distribution of a over the 
SNR show spatial variations which steepen towards the inner western filamentary shell. Polarized emission is prominent on the 
western filamentary shell region. The RM there is estimated to be about 30±77n rad m"-, where the n=l solution is preferred, and 
the magnetic field has a strength of about 9.5 /jG. From the HI channel maps, further evidence is provided for an interaction between 
the Monoceros SNR and the Rosette Nebula. We identify partial neutral hydrogen shell structures in the northwest region at LSR 
velocities of +15 km s"' circumscribing the continuum emission. The HI shell has swept up a mass of about 4000 Mq for a distance 
of 1.6 kpc. The western HI shell, well associated with the dust emission, is found to lie outside of the radio shell. We suggest that the 
Monoceros SNR is evolving within a cavity blown-out by the progenitor, and has triggered part of the star formation in the Rosette 
Nebula. 

Conclusions. The Monoceros SNR is interacting with the ambient interstellar medium with ultra-high energy emission detected. 
Its interaction with the Rosette Nebula is further supported by new evidence from HI data, which will help the investigation of the 
emission mechanism of the high energy emission. 

Key words. ISM: supernova remnants - Polarization - Radio continuum: general - Methods: observational 
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1. Introduction 

The Monoceros Nebula (SNR G205.5+0.5) is an extended fil- 
amentary object lying in the anti-center complex region. The 
bright Rosette Nebula (SH 2-275) and HII region SH 2-273 are 
adjacent to the edge of the southern and northeastern shell of the 
SN R. It w a s first recognized as a possible supernova remnant 
bv iDaviesI (Il963h from radio observations at 237 MHz. More 
detailed observations presented bv lHoldeiil (1 1 9681) confirmed this 
property by a non-thermal radio spectrum. 

The SNR has be en stud ie d at l ow and intermediate 
radio frequencies by 'Holdten| (Il968l) : iMilne & Hilll d 19691) 
and Dickel & Denoy er ( 1975 |). Observations a t frequencies near 



2700 MH z were made by Dav et alJ (Il972h : 'Milne & Dickel' 



1974); 



, .Velusamv & Kundu ( 1974b and Graham et al.. 

19821) . iGraham et al.1 (1 19821) found an averag spectral in- 
dex a = -0.47 + 0.06 from 111-2700 MHz radio data. 
They also presented the spatial spectral index distribution at 
178-2700 MHz and 611-2700 MHz, which was found to be 
steeper in the inner region. They also identified an Hll region 
(G20 6. 35+1.35) in the eastern periphery from its thermal spec- 
trum. 'G ao et al" I (12011') showed a new polarization observation 
at 4800 MHz. Combined this with Effelsberg 21 cm and 1 1 cm 
data, they obtained a spectral index of a - -0.43 + 0.12 for 
the western shell. Early 1 1 cm polarization observations were 
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made by iMilne & Dickell d 19741) . but no significant polarized 
emission was detected. The 6 cm polarization map showed that 
the magnetic field is closely aligned with the filamentary shell. 

The distance to the Monoceros SNR has been estimated 
to range from around 800 pc to about 1.6 kpc. lOdegardI 
(1986) summarized all the available evidence and argued that 
the Monoceros SNR is at the same distance as the Monoceros 
OB2 association of about 1.6 kpc. This is further supported 
by the decimeter observations of the absorption of non- 
thermal emission from the supernova remnant by the Rosette 
Nebula (Odegard 1986). The diameter of the SNR is then about 
106 pc (3.8°), an d the age is about 1 .5 x 10^ yr A study of 
Einstein IPC data (iLeahv et al.ll985Lll986h shows diffuse X-ray 
emission from the Monoceros SNR shock front, coiTesponding 
to the optical filaments. By fitting the data, they concluded that 
the SNR was expanding in a low-density (0. 003 cm'-^) medium , 
and ha s not y et entered the radiative phase (iLeahv et alj|1986l) . 
lRiddlel(ll972h identified a neutral hydrogen shell located outside 
of the Monoceros optical filament with low resolution HI data. 

The Monoceros SNR is suggested to be interacting with 
the Ro sette Nebulae. I n the SNR ridges overlapping Rosette 
region, [Fountain et al. I (fT97l found broadened line widths of 
H109ff and much enhanced Ha emission than average in 
the Rosette Nebula. EGRET detected extended y-ray emis- 
sion (3EG J0634+0521) from the region associated with the 
Monoceros/Rosette Nebula in the energy range from 100 MeV 
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up to 10 GeV (iJaffeet al.lll9"97h . This emission was interpreted 
as y-rays from the decay of tt^'s produced by the interaction of 
shock accelerated protons with the ambient matter. HESS found 
a source HESS J0632+057 (1 00 GeV - 10 TeV) located at the 
rim of the western SNR shell (iFiasson et al ]l2008 l). which seems 
to be associated with 3EG J0634+0521 and an associated molec- 
ular cloud. An observation carried out with the BeppoSAX satel- 
lite (Kaaret et al. 1999) discovered a hard spectrum X-ray point 
source (SAX J0635-I-0533) within the 95% probability circle of 
the EG RET detection, which was later iden tified as a binary 
pulsar dCusumano et al. I l2000t) . [Torres et al.l (2003) postulated 
that both the binary pulsar SAX J0635+0533 and the emission 
from the interaction region between the expanding shell and the 
Rosette Nebula contribute to the EGRET y-ray emission. 

Although there is akeady a considerable body of radio ob- 
servations on the Monoceros SNR, we try to better delineate the 
extent of the SNR and the spatial spectral index variations using 
new radio data with higher resolution. Also, a new high sensitiv- 
ity, high resolution HI data cube enables us to identify HI fea- 
tures morphologically associated with the SNR and to discuss 
various evolutionary scenarios. In Sect.2, we describe the 11, 
21, and 6 cm radio continuum observations, as well as the HI 
data of the Monoceros SNR. The determination of the flux den- 
sities and the spectral index distribution are presented in Sect. 3. 
The polarization properties and the magnetic fields in the north- 
east filamentary shell are estimated in Sect.4. The associated HI 
shell analysis is presented in Sect.5. Finally, we give a discussion 
and summary in Sects. 6 and 7, respectively. 



2. Data 

2.1. Radio data 

The 21 cm and 1 1 cm radio continuum data presen ted here are 
taken from the EfFelsberg Galactic plane survey dReich et al.l 
[1990; Furst et al. 1990). The angular resolution of these sur- 
veys are 9.4' and 4.4', respectively. The 6 cm continuum and 
polarization data are extracted fro m the Sino-Germ an polariza- 
tion survey of the Galactic plane (iGao et al ] l2010l) made with 
the Urumqi 25-m telescope. The angular resolution is 9.5' at 
4.8 GHz. Detailed system and observation parameters and the 
data reduction procedure are described in these related papers. 

Figure [U shows the Effelsberg 21 cm and 11 cm total in- 
tensity images of the Monoceros SNR. The Urumqi 6 cm con- 
tinuum emission with corresponding polarization intensity (E- 
vectors) overlaid are presented in Fig.|2] Filamentary shell struc- 
tures are clearly seen in the northern, western and southern re- 
gions with diffuse emission in the eastern portion. The HII region 
0641-H06 (/=206?35, /7=1°35) was identified by Graham etaD 
(Il982h . There is also a small emission ridge (/ = 2QT.6,b - 
-1?3) projecting from the Rosette Nebula towards the east. In 
Sects. 3 and 5, we propose this structure to be part of the 
Monoceros southern shell. 

2.2. HI data 

The HI data sets come from the recently released Galactic 
Arecibo L- Band Feed Array HI survey (GALFA-HI), described 
in detail bv P eek et al.l(l201 \ \. The angular resolution of the sur- 
vey is about 4'. It covers a wide velocity range from -700 km s"' 
to -1-700 km s"' with a velocity resolution of 0.18 km s"'. Typical 
noise levels are 80 mK in an integrated 1 km s ' channel. 
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Fig. 1. The Effelsberg 21 cm and 11 cm continuum images of 
the Monoceros SNR, with angular resolutions of 9.4' and 4.4' 
respectively. Contours at 21 cm start from 100 mK Tb and in- 
crease by 2"x80 mK Tb (n=0,l,2,3,- ■ •)• Contours at 11 cm are 
drawn according to 2"x30 mK Tb, with n=0,l,2,3,- ■ ■. 



3. Spectral index analysis 

3. 1 . Subtracting the thermai radio continuum emission 

Thermal emission is well traced by the infrared (IR) emis- 
sion from the dust. We present the IRIS 60 jum im- 
age of the Monoceros complex region in Figure [3] The 
IR observations are the high-resolution reprocessed IRAS 
image produce d at the Infrared Processing and Analysis 
Center (IPAC) dMiviUe-Deschenes & Lagachd l2005l) . Besides 
the strong thermal emission from the Rosette Nebula, there is 
a thermal source (/=206?25, b^(f.75) within the SNR region and 
HII region 0641+06 within the shefl. 

We tried to separate and subtract the aforementioned ther- 
mal sources using the correlation between the radio bright- 
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Fig. 2. The Urumqi 6 cm total intensity map of the Monoceros 
SNR at an angular resolution of 9.5'. Polarization vectors in the 
B-field direction (E+90°) are overlaid. The length of the bars is 
proportional to PI, with a polarized intensity of 1 mK Tb corre- 
sponding to a bar-length of 2'. Total intensity contours are drawn 
according to 2"xl0 mK Tb, with n=0,l,2,3,- ■ -. 
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Fig. 3. IRIS 60 fim emission of the same area as in the preceding 
radio continuum images. Resolution is 2'. Contours of 6 cm total 
intensity are overlaid having the same levels as for Fig. |2] 



ness temperature and its equivalent dust fl ux at 60 /im 
(e.g. lBroadbentet"an (Il989l) : iFiirst et al.1 (I1987I) '). This correla- 
tion generally exists for diffuse thermal emission and has been 
used to remove the thermal background emission from SNRs. 

For diffuse HII regions, Broadbent et al. (1989') obtained the 
relation Tb,\\ = (6.4 + 1 .7) x /eo, where Tb.i i is in mK and /eo in 
MJy sr"' . Assuming a temperature spectral index /3 = -2. 10, the 
corresponding relation at 21 cm becomes T^.n = 24.6/60. We 



made a TT-plot between 21 cm and 6 cm for the thermal source 
region. A temperature spectral index /3 = -2.10 + 0.22 confirms 
the thermal nature of the emission. We convolved the 21 cm and 
60 fim images to the same 10' resolution, and plotted the bright- 
ness at 60 fim (/eo) against that at 21 cm (rB,2i). The resulting 
correlation is shown in Fig. |4] from which we find a slope of 
23.6 mK MJy"' sr, comparable with the theoretical value. We 
were then able to estimate the theoretical thermal emission and 
simply subtract it from the raw images. The thermal emission at 
1 1 cm and 6 cm are simply interpolated using a spectral index of 
/?= -2.10. 

We checked the similar relation for the Rosette Nebula. Both 
the Effelsberg 21 cm total intensity data and the IRIS 60 yum 
data were convolved to 10' and scaled to units of Jy/beam. The 
ratio Ts,2\/ho shows large scatter in the periphery, with an av- 
erage value of about 90 mK MJy"' sr in the central region. In 
regions where the two objects overlap each other, it is difficult 
to separate the nonthermal emission from the SNR and the ther- 
mal emission from the Rosette Nebula. Therefore we directly 
exclude the Rosette Nebula emission by setting a boundary in 
between to obtain the emission from the SNR alone. 
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Fig. 4. Plot of the brightness temperature T^.n against the bright- 
ness /(,(). The best-fit straight line corresponds to a slope of 
23.6 mK MJy-' sr. 



3.2. Integrated flux densities and TT-plots 

In order to estimate integrated flux densities and derive an ac- 
curate spectral index for the remnant, we removed five bright 
point- like sources from the 21,11 and 6 cm total intensity maps 
within the area of the Monoceros SNR by Gaussian fitting. 

We use d the method of tem perature-versus-temperature plots 
(TT-plots) dTurtle et alj|1962 ) to adjust the base-levels for the 
entire SNR area at two wavelengths. Both maps were convolved 
to a common angular resolution of 10'. Corresponding to the TT- 
plot results for the entire shell region, we found a constant offset 
of 3 mK Tb at 6 cm, 55 mK Tb at 1 1 cm and about 35 mK Tb at 
21 cm. 

The integrated flux densities are obtained by setting poly- 
gons just outside the periphery of the SNR and integrating the 
enclosed emission. The contribution of the Rosette Nebula was 
excluded by setting the highest brightness temperature in the 
shell as its boundary. From variations outside the SNR, we 
estimated a remaining uncertainty in setting the base-level of 
9.9 mK Tb at 21 cm, 7.9 mK Tb at 11 cm, and 3.5 mK Tb 
at 6 cm. We obtained a flux density of 120.8+9.9 Jy for the 
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Fig. 5. Spectral index from integrated flux densities of the 
Monoceros SNR. 



Monoceros SNR at 21 cm, 89.7+7.9 Jy at 1 1 cm and 73.4+3.5 Jy 
at 6 cm. The linear fit to these integrated flux densities for the 
Monoceros SNR is shown in Fig.|5] yielding a spectral index of 
a = -0.41+0.16. 

The TT-plot method is unaffected by the uncertainty of the 
background levels in the maps. We also applied it to an investiga- 
tion of the spectral index for different sections of the Monoceros 
SNR shell. The temperature spectral index p {J3 - a - 2) found 
from fitting the slope of 2 1 cm and 6 cm pairs of different shell 
regions are shown in Fig. |6l The error of y6 for the diffuse east- 
ern shell is relatively large, probably due to its low brightness 
and confusion with weak unresolved background sources which 
could not be subtracted. We also investigated the southern shell 
region which is mixed with the Rosette Nebula. Beside the old 
spherical shell region, there is the new shell-like feature protrud- 
ing out towards the east (see Section 2.1). The index of the new 
shell-like feature is about a - -0.53 + 0.44. The large uncer- 
tainty is due to the weak emission, plus the residuals after a 
point-source has been subtracted. Ho wever, the spe c tral in dex 
distribution in this region in Fig. 4 of iGraham et al.l ( |1982|) . as 
well as Fig. 7 in Sect. 3.3, are compatible with the feature being 
nonthermal. It seems that this feature is likely to be a part of the 
Monoceros SNR. 

3.3. Spectral index map 

We calculated the spectral index for each pixel of the Monoceros 
SNR from the zero-level corrected images at the three frequen- 
cies, by linearly fitting for intensities versus frequencies using a 
logarithmic scale. In order to achieve reasonable spectral indices 
without significant influence from noise or local distortions, a 
lower intensity limit of 180 mK Tb, 30 mK Tb and 6 mK Tb 
are set for the 21, 11 and 6 cm maps, respectively. The regions 
outside the remnant are blanked out. We display the tempera- 
ture spectral index map from 21 cm, 1 1 cm and 6 cm in Fig. Q 
Possible remaining variations of the base-levels at 1 1 cm or 6 cm 
yield an average systematic uncertainty of the spectral indices of 
Aa ~ 0.2 with errors being larger in regions of weak emission. 

The Rosette Nebula shows a flat brightness spectral index, as 
expected for optically thin thermal emission from an HII region. 
The spectral index of the SNR G206.9+2.3 shows fairly uniform 
distribution and has an averaged index of yS = -2.46, similar 
to the value found for its integrated flux densities by Gao_e t_al., 
( 1201 lb . Spectral index values are generally about -0.4 in the 
Monoceros SNR shell region. The new shell-like feature in the 



southern shell has an index of about -0.5. Little significant spec- 
tral variation is traced in the overlap region with the Rosette 
Nebula. The spectral index gets steeper in the inner region so 
that the inner western filament at / = 205?5, b - 1°0 has a steeper 
spectral index of about -0.6. The value is generally consistent 
with that in the l ower frequency spec tral index map between 178 
and 2700 MHz (lOraham et al.lll982h . 




I 1 \ \ 1 1 
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Fig. 7. Temperature spectral index map of the Monoceros SNR 
from 21/1 1/6 cm maps. Regions where the random error exceeds 
0.06 have been blanked out. Contours of 6 cm total intensity with 
the same levels as Fig.|2]are overlaid. 



iGraham et al.l ( Il982h excluded the possibility that the differ- 
ence comes from a mixture of thermal emission for the more 
northerly rim of NGC 2264 using a positional argument and by 
comparison with the lower frequency index. The new higher fre- 
quency 6 cm data further excludes this possibility. The relativis- 
tic electrons in the inner western ri dge probably have a d iffer- 
ent energy spectrum as proposed by IGraham et al.l (Il982h . It is 
now generally accepted that relativistic electrons responsible for 
the synchrotron emis sion are produced by the Fermi accelera- 
tion mechanism ( Belilll978afbl) . This process predicts a spectral 
index given by -/3 = 2 + (M^ + 3)/[2(M^ - 1)], where M is the 
Mach number of the shock. For a very strong shock, M is close 
to 3, corresponding to an average value of /3 - -2.75. A weaker 
shock produces a steeper spectrum. The shock in the inner west- 
ern filament is possibly weaker with a decreased compression ra- 
tio and results in a steeper spectrum. In addition, a steeper spec- 
trum could also have originated from high-energy electrons due 
to synchrotron aging, when propagating in a weak magnetic field 
region. However, the strength of the magnetic field in both west- 
ern ridges appears to be similar from the vectors in the 6 cm map, 
both from the perpendicular component and the amounts of 
Faraday rotation. No indication of a spectral break is found in the 
frequency range from 178 MHz to 4.8 GHz, which means that 
any spectral break due to synchrotron aging should be at higher 
frequencies. Higher frequency data are needed to distinguish be- 
tween these two explanations. 
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Fig. 6. TT-plots between Urumqi 6 cm data and the EfFelsberg 21 cm data (as examples) for the eastern shell, northwest filamentary 
shell and newly identified southern branch shell of the Monoceros SNR. 



4. Polarization analysis 

4.1. The linear polarization properties 

The 6 cm polarization map is shown in Fig. |2] Prominent polar- 
ized emission is mainly displayed towards the western shell of 
the SNR. The observed magnetic vectors are nearly distributed 
along the ridges of the two branches. This indicates a well or- 
dered magnetic field in th e shell, as expected for an evolved 
SNR (iFiirst & Reichll2004l) . The polarization percentage of the 
inner filament is about 20% at 6 cm, and less at about 10% in the 
outer filament. Except for the enhanced polarization along the 
SNR's western shell, polarized emission in other shell regions 
is very weak. The inner polarization patch might originate from 
polarized diff'use interstellar emission possibly mixed with polar- 
ized emission from the SNR of about similar strength. Strongly 
polarized emission is also detected from the SNR G206.9-H2.3. 

The large-scale emission component is missing in the 6 cm 
polarization map, due to the zero-level setting during observa- 
tion and data analysis. This is known to limit the interpreta- 
tion of polarization structures cau sed by Faraday effects in the 
interstellar medium dReichl 12006 ) except in strongly polarized 
regions. We checked the polarization map with large-scale po- 
larized emission added, which is extrapolated from the WMAP 
K-band data usi ng an appropriate spectral index (see Fig. 6 in 
iGao et all (l2010h V Strongly polarized emission from the diffuse 
interstellar medium appears near the Galactic plane. The polar- 
ized signal from the western shell of the Monoceros SNR is suf- 
ficiently strong that it remains almost unchanged by adding the 
large-scale emission component. 

We extracted the 21 cm polarized intensity map of the 
same region (Fig. [8]) from the polarization survey made with 
the DRAO 26-m telescope ( Wolleben et al. 2006). This map in- 
cludes the large-scale polarized emission component at an abso- 
lute zero level with an angular resolution of 36' . We found a large 
diffuse interstellar polarization patch covering the Monoceros 
SNR region, with enhanced polarized emission coinciding with 
the western filamentary shell region. The polarization angles in 
the western shell lie nearly perpendicular to the shell. Polarized 
emission from a larger distance might has been depolarized to- 
wards the Rosette Nebula, where weakly polarized emission is 
detected at 21 cm. 

Compared to the western shell, the low polarization percent- 
age of the eastern shell indicates a more irregular magnetic field. 
The background polarization there is depolarized by small-scale 
magneto-ionic fluctuations in the shell and/or polarized emis- 
sion from the SNR shell with a different orientation. In the 
southern shell where interaction has occurred, as discussed in 
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Fig. 8. DRAO 21 cm polarized intensity map of the Monoceros 
SNR at an angular resolution of 36'. Polarization vectors in the 
B-field direction (E+90°) are overlaid. The length of the vec- 
tors is proportional to PL A polarized intensity of 1 mK Tb 
corresponds to a length of 0.05'. Contours show total intensi- 
ties starting at 200 mK Tb and increasing by 2" x60 mK Tb 
(n=0,l,2,3,- ■ ■)• The positions of RM sources in iTavlor et all 
(2009) are indicated by crosses, with lengths proportional to the 
RM values. 



Sect. 5.3, the high thermal electrons densities in the Rosette 
Nebula (Odegard 1986) and their fluctuations could be suffi- 
ciently strong to cause a heavy depolarization of the background 
polarization. 

4.2. Rotation measures of the western shell 

Since the polarized emission from the western shell is strong and 
clearly stands out against the background in the 21 cm polariza- 
tion map, we were able to remove the diffuse Galactic polarized 
emission component and calculate the rotation measure (RM) 
distribution using the polarization angles from the 21 and 6 cm 
polarization maps. We subtracted an average value of 81 mK Tb 
from the 21 cm U map and "twisted" the Q map setting by values 
at the four corners; (from upper left to lower left these values are 
8, 80, 4, 8 mK Tb). The 6 cm U and Q data were smoothed to an 
angular resolution of 36' and we re-derived a map of polariza- 
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tion angles. We linearly fit the polarization angles at 21 cm and 
6 cm versus the square of the wavelength for each pixel from all 
maps. To gain a high signal to noise ratio, pixels with brightness 
temperatures less than 3xcr/)/ were not included. 

In principle, we require polarization maps for at least three 
frequencies to obtain unambiguous RMs. RMs calculated at 
21 cm and 6 cm have an ambiguity of +n x 77 rad nT^ (n=0,+ l, 
+2,...). The average minimum RM (n=0) is about +30 rad m"^ 
in the western filamentary region. Larger RMs are also possi- 
ble in this case. We checked the RMs of extragalactic sources 
towards this region from the catalog by Taylor et al. (2009). 
There is an extragalactic source (/=204°92, b=Q°.86) with an 
RM of 202.1+6.6 rad m"^ shinning through this region, and 
another source (/=205°02, b = -0°96) lies in the southeast 
shell with an RM of 152.7+5.4 rad m"^. Two sources within 
the Monoceros SNR region have RMs of 47.1 ±1.8 rad m"^ and 
58.3+2.2 rad nT^ respectively (see Fig.[8]l. Though the RMs of 
the extragalactic sources represent RMs integrated through the 
total interstellar medium, the RM values of about 109 rad nT^ 
(n=l) and 184 rad nT^ {n-2) are more likely from the western 
shell, considering the distance of the Monoceros SNR of 1 .6 kpc. 
In Fig. [1] we noticed that the polarization B vectors at 6 cm in the 
western shell as shown are distributed almost tangential to the 
shell. If we accept deviations of about 20° from an exactly tan- 
gential orientation, the RM is estimated to be about 96 rad m"^, 
near to the RM value for n = 1 of 109 rad nT^. Detailed inves- 
tigation of the unambiguous RMs requires polarization observa- 
tions at other wavelengths. 

4.3. Strength of the magnetic field in the western shell 

We can estimate the magnetic field strength in the western shell 
by assuming energy equipartition between the magnetic field 
and electrons and protons in the SNR-shell. We used the re- 
vised equipartition formula by Beck & Krause (2005), which 
uses the number density ratio K instead of the total energies of 
cosmic ray protons and electrons in classical equipartition es- 
timates. Assuming the synchrotron emission is concentrated in 
the SNR filamentary shell, the radio emissivity Ly/V averaged 
over the source's volume is replaced by the surface brightness 
over the line-of-sight path length through the emitting region 
ly/l. Therefore: 

_ 4n(-2a+l){K, + l)IyEj;^''(v/2cir\i/(-a^3) 



where, ly is the synchrotron intensity at frequency v and a the 
synchrotron spectral index of -0.4. Ko is a constant ratio of num- 
ber densities of protons and electrons in the energy range, which 
is assumed to be 1. Ep is the pr oton rest energy, ci a nd C2 are 
constants (see the Appendix of Bec k & Krausd (l2005h '). We in- 
tegrated the flux density at 6 cm for the western shell section 
centered at I = 204?6, b = Q°5 {{6l,6b) = 2?3 x 0^8), and ob- 
tained a value of 1 1 Jy. The radio intensity of the western shell 
is about ly - 0.2 x 10^'^erg s^'cm^^Hz^'sr"'. We estimate the 
thickness of the western filament from high angular resolution 
11 cm data to be ~ 13', about 10% of the radius of the SNR. 
The thickness of the filament is about 6 pc, and the total path 
through the interacting zone is / ~ 49 pc at a distance of 1.6 kpc. 
The estimated equipartition magnetic field Bgq is about 9.5 yuG, 
corresponding to a factor of ~ 3-4 enh ancement of the mean 
ISM field of B,o, 2-3//G (iHan et alJl2006h . 



5. Neutral hydrogen morphology 

5.1. HI associated with the Monoceros SNR 

The GALFA-HI survey covers a velocity range from -700 km 
s~' to +700 km s"' with a velocity resolution of 0.18 km s"'. 
iRaimondl (1 19661) and|Riddle (1972) have shown that a neutral 
hydrogen shell is located outside of the Monoceros optical fila- 
ment with an HI cavity inside. The GALFA-HI image has much 
better resolution and sensitivity and reveals many detailed struc- 
tures. 

The 0.8-1.6 kpc distance corresponds to a velocity range 
o f +8 - 15 km s ~' according to the Galactic rotation model 
( (lFichetal.lll989h : Rq = 8.5 kpc and 0o = 220 km s '). We 
have inspected the entire data cube of this complex region from 
the GALFA-HI survey in the velocity range between km s"' 
and +30 km s"', looking for structures morphologically related 
to the SNR. 

Fig. |9] displays the HI channel maps after averaging over 
11 consecutive channels (~2 km s"') within the velocity inter- 
val from +4 to +25 km s '. The central velocity of each im- 
age is indicated at the top. Contours of 6 cm total intensity 
(smoothed to 10') are overlaid to indicate the Monoceros com- 
plex region. Towards this region, neutral hydrogen filamentary 
structures caused by stellar winds from the OB associations of 
Mon OBI and Mon OB2 are visible across the entire region. The 
large-scale distribution of neutral hydrogen gas shows a den- 
sity gradient towards the Rosette Molecular Complex (RMC). 
The HI shell related to the RMC lies between +5 km s ' and 
+25 km s"', with a lower emission towards the southern shell 
of the Monoceros SNR, which is consistent with the result 
in lKuchar & Banial(ll993h . 

We have identified a partial neutral hydrogen shell surround- 
ing the western SNR, which we believe is associated with the 
SNR. This shell structure is distributed over velocities between 
+7 km s"' and +22 km s"' with a thickness of ~13'. It lies 
slightly outside the spherical outline of the remnant, in partic- 
ular in the southwest direction using Galactic coordinates, and 
seems to have expanded ahead of the radio shell. This is similar 
situation to the case of another likely o ld SNR, the North Polar 
Spur (Loop I) (Berkhu iisen et al.llT970l) . However, the HI shell 
is well associated with the 60 yum emission (Fig. |3]l, suggesting 
that this shell could have been blown-out by strong stellar winds 
from the progenitor of the SNR before the explosion. The 21 cm 
line emission is weak in the eastern SNR sector It seems that the 
SNR blast wave has expanded into a medium of lower density in 
this direction. In the southern SNR region, the HI shell features 
are mixed with those of the Rosette Nebula. However, a small 
section of possible shell structure is seen in the HI emission for 
the velocity range +10 to 14 km s"', projecting from the east 
of the Rosette Nebula. This is prominent, being 0°5 long with 
a thickness of 10' in the +12 km s ' channel map, and is well 
associated with the new shell-like radio continuum feature in the 
southern shell. 

In Fig.[TO]we present the result for the integrated HI column 
densities over the velocity range +4 to +25 km s Weak HI 
emission segments associated with the eastern radio shell of the 
Monoceros SNR emerge after the large-scale diffuse emission 
has been subtracte d using t he "un sharp-masking" procedure de- 
scribed by So fue & ReichI (Il979l) . A cavity with an angular di- 
ameter of about 3.8° is visible coinciding with the Monoceros 
SNR. 

Assuming the HI emission to be optically thin, we calculate 
the average column density for the northwest shell to be about 
5.5 X 10^" cm"^. If a physical association exists, with a distance 
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Right Ascension fjaOOO) Rif^ht Ascension (J2000) Risht Ascension fjSOOO^ 

Fig. 9. HI channel-maps in the area of the Monoceros SNR for the velocity range from +4 km s ' to +25 km s \ averaged over 1 1 
consecutive channels, yielding a velocity resolution of ~2 km s"'. The central LSR velocities are indicated above each panel. Total 
intensities of 6 cm continuum emission smoothed to 10' are overlaid as contours. 
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Fig. 10. Integrated Hl-intensities for the velocity range from 
+4 km s"' to +25 km s"' with contours of 6 cm continuum emis- 
sion overlaid, having levels of 10, 30, 120, 640 mK. 
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Fig. 11. Integrated CO-intensities for the velocity range from 
km s"' to +30 km s"' with 6 cm continuum contours overlaid, 
having levels of 10, 30, 120, 640 mK. 



has a ring shape around the nebula with a gap towards the 
Monoceros SNR. Based on the associated CO morphology and 
the fact that the Monoceros SNR has a similar velocity to 
that of the Rosette Molecular Complex, we suggest that the 
Monoceros SNR is interacting with the Rosette Complex, and 
that the progenitor of the Monoceros SNR comes from the 
Rosette Complex. 



of 1 .6 kpc, the thickness of the half-shell of 1 3' with a 1 .9° radius 
corresponds to a swept-up mass (wihA^CHI)) of about 4140 M©, 
assuming the southern shell has a similar column density to the 
western shell. The average density in the Hl-shell is then about 
n ~ {NIL) = 3.6 cm"^ with a line-of-sight depth within the shell 
of 49 pc. If this mass of neutral gas was originally uniformly dis- 
tributed within a sphere of radius 1 .9°, the pre-explosion ambient 
density would have been 0.27 cm""*, consistent with the ambient 
density of 0.34 cm"^ estimated by Odegard ( 1986), according to 
the Chevalier formula (Che valieitil974 ). The evacuation by the 
action of stellar winds within the Mon OB2 association and the 
progenitor may cause the low density environment. 

5.2. CO associated with fhe Monoceros SNR 

We inspected the CO images of this region obtained 
by [Oliver et akj (Il996h . There are weak CO clouds in the direc- 
tion of the enhanced radio continuum emission of the western 
shell with velocities ranging from -3.1 km s"' to +9.4 km s"', 
and enhanced CO emission is found towards the southern shell 
with velocities ranging from + 10 to +30 km s"' . The Cone cloud 
adjacent to the northwest of the Monoceros SNR has an average 
velocity of +7 km s"'. This cloud is the progenitor of the young 
open cluster NGC 226 4 and belongs to the Mon OB 1 associa- 
tion (I01iveretal.lll996h . 

In Fig. [TT| we sh ow the integrated '^CO emission 
from lOliver et aP (1 19961) in the velocity interval < v < 
+30 km s"' with contours of the 6 cm continuum emission over- 
laid. The molecular cloud associated with the Rosette Complex 



5.3. Interaction witli tlie Rosette Nebula 

The high resolution HI channel maps provide evidence that the 
Monoceros SNR is interacting with the Rosette Nebula. In the 
integrated HI emission map of Fig. [10] the Rosette Nebula ap- 
pears to have broken into the Monoceros SNR area outlined by 
its outer HI shell, the HI shell associated with the Rosette Nebula 
presents a ring shape with a minimum in the direction of the 
Monoceros SNR. In the channel maps of Fig.|9] a small section 
of possible HI shell structure is seen at v — hl2 km s"' project- 
ing from the east of the Rosette Nebula, coinciding closely with 
the possible new radio shell. This feature might have a relation- 
ship with the interaction between the Monoceros SNR and the 
Rosette Nebula, and is possibly caused by the shock of the SNR 
or stellar winds in the Rosette Nebula. 

It has been previously speculated that the formation of 
the cluster NGC 2244 was triggered by an encounter betwee n 
the Monoceros SNR and the RMC dTownslev et alj l2003h . 
Morphologically, it seems that the progenitor of the SNR ex- 
ploded first and its shock formed a well-shaped shell. The ex- 
pansion of the Rosette Nebula ionized by the cluster NGC 2244 
collided with the SNR shell and kept expanding. However, the 
SNR's age has been estimated to be 0.3-5x10^ yr (^ Welsh et aP 
2001). Th e main-seque nce turnoff age of the cluster is about 
1.9 Myr (lPark& Sungll2002h . which suggests that the cluster 
formed earlier than any possible encounter with the SNR. The 
kinetic energy of t he expansion of the R osette Nebula was about 
~ 3.8 X 10"*^ ergs ( Kucliar & Banial[T993l) . which is smaller than 
the typical SNR explosion energy of 10^' ergs. Although there is 
little nonthermal radio emission detected in the Rosette Nebula, 
it is still possible that the shock of the Monoceros SNR has bro- 
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Fig. 12. Upper panel: Ha emission overlaid with 1 1 cm continuum contours having the same levels as for Fig. [T] Lower panel: 
ROSAT X-ray emission smoothed to 10' angular resolution overlaid with 6 cm total intensity contours haveing the same levels as 
for Fig. m 



ken into the Rosette Nebula at some places, possibly the new 
radio shell region, and is triggering new star formation in the 
cluster. 

Despite the fact that no detection of maser emission has bee n 
reported in the radio survey at 1720.5 MHz bv lFrail et an(ll996h . 
the Monoceros SNR has been interpreted as interacting with the 
atomic and molecular components of the Rosette Complex envi- 
ron ment based on enhance d line widths of HlOQa- and Ha emis- 
sion ([Fountain et al.|[T979l). and similar velocity of Ha emission 
in the contact region (lRevnoldsl ll983). EGRET and HESS have 
detected high energy y-ray emission from the overlap region be- 
tween the Monoceros SNR and the Rosette Nebula, interpreted 
in terms of tt" decay resulting from collision s between SNR- 
accelerated partic les and the molecular cloud dJaffe et al.lll997t 
[Torres et al.l2003h . For the first time, the present HI observations 
provide direct morphological evidence of their interaction. The 
Monoceros SNR has swept up an HI mass of about 4000 M©, 
with a possible break in the shell at -1-12 km s"' towards the 
Rosette Nebula. A number of Galactic SNRs in ph ysical con- 
tact wi th molecular clouds have been summarized bv lJiang et al.l 
dMlO'). A similar case is the SNR HB3, which is interacting with 
the adjacent HII region W3 and the associated star-formation re- 
gion dRoutledge et al.lll99lh . 



6. Comparison with other multi-bands images 

We compare multi-wavelength observations to discuss the phys- 
ical properties and the environment of the Monoceros SNR. 

The Monoceros SNR shows diffuse Ha emission associated 
with most of the SNR shell region, and outside of the shell to the 
southwest shell (Fig. fT2l i. The diffuse emission in the area east 
of the SNR represents a low density environment, which results 
in a less compressed magnetic field and thus reduced efficiency 
of particle acceleration. The radi o filaments are wel l correlated 
with the optical filaments ([Nil]) (iDavies et al .119781) . except for 
the inner western filament. This correlation is not seen in the 



optical image, probably due to lower temperatures/densities in 
the radio filaments. 

The Monoceros IRIS 60 fim dust emission with 6 cm 
total intensity contours superposed is shown in Fig. [3] The 
shell-associated dust emission stretching from the eastern spur 
of the Rosette Nebula is enhanced with high ISM density. 
The dust emission envelope in the western region is out- 
side the radio shell, and is associated with the diffuse optical 
shell (Davies et al. 1978). It is also well circumscribed by the 
HI shell, suggesting that the dust there has been ionized by a 
pre-shock caused by stellar winds from the progenitors. 

The Monoceros image (0 .56-1.21 keV) f rom the ROSAT soft 
X-ray all-sky survey ( Snow den et al.|[T997l) is shown in Fig. [12] 
with contours of 6 cm total intensity overlaid. Diffuse X-ray 
emissions is located on the rim of the remnant and coincides 
with the densest region of optical filaments, as expected for a 
young SNR in the Sedov phase. In the eastern SNR shell, the 
X-ray emission is located outside the boundary of the radio con- 
tours which are assumed to trace the outer SNR shock, indicating 
mo re diffuse emission in the weak eastern shell. 

iLeahv et af] (ll985l[T986l) fitted the weak X-ray emission de- 
tected by the Einstein X-ray satellite, and concluded that the ini- 
tial density of the medium in which the Monoceros supernova 
occurred was 0.003 cm"^, and the expansion of the SNR pro- 
ceeds in a non-homogeneous multi-component medium. The HI 
environment confirms this and shows a gradient of density. The 
Monoceros SNR has probably expanded in a pre-existing cav- 
ity generated by the progenitor. In this case, the cavity would 
have had to be extremely large, about 105 pc, which may have 
required more than one star to make. The shock has recently en- 
countered the interstellar material, which turns it into a bright 
emitter at radio wavelengths. 

7. Summary 

After subtraction of the inner thermal emission and compact 
extragalactic sources, we obtained integrated flux densities of 
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120.8±9.9 Jy at 21 cm, 89.7±7.9 Jy at 11 cm and 73.4+3.5 Jy 
at 6 cm for the Monoceros SNR. The spectral index from these 
three wavelengths is about a = -0.41 + 0.16, which was con- 
firmed by TT-plots. A new southern shell branch is found, likely 
belonging to the Monoceros SNR, based on its spectral index of 
~ -0.5. The distribution of spectral index over the SNR shows 
some variations, and steepens towards the western inner filamen- 
tary region, corresponding to a weaker shock compression. 

Strong polarized emission is observed from the western fila- 
ment of the Monoceros SNR at 6 cm and 21 cm, indicating the 
presence of a strong regular magnetic field component. The RM 
values there are ambiguous, and are typically 30±77n rad m"^, 
with n = 1 being the favored value. The magnetic field in the 
western shell region is estimated to be 9.5 yuG. 

The case for interaction between the Monoceros SNR and 
the Rosette Nebula is strengthened by the HI data. We iden- 
tify partial neutral hydrogen shell structures at LSR velocities 
of +15 km s"' circumscribing the continuum emission. The HI 
shell has swept up a mass of about 4000 M© for a distance of 
1.6 kpc. The Monoceros SNR has probably triggered part of the 
star formation in the Rosette Nebula. The western HI shell is 
found to be well correlated with the dust emission outside of 
the radio shell, indicating that the Monoceros SNR is evolving 
within a cavity blown-out by the progenitor, and that the SNR is 
still in the Sedov phase. 
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